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Abstract 

In recent years, much study has been made by applying the non- extensive statistics 
(NES) to various non-extensive systems where the entropy and/or energy are not neces- 
sarily proportional to the number of their constituent subsystems. The non-extensivity 
may be realized in many systems such as physical, chemical and biological ones, and also 
in small-scale nanosystems. 

After briefly reviewing the recent development in nanomagnetism and the NES, I 
have discussed, in this article, NES calculations of thermodynamical properties of a nan- 
ocluster containing noninteracting M dimers. With bearing in mind a transition-metal 
nanocluster, each of the dimers is assume to be described by the two-site Hubbard model 
(a Hubbard dimer). The temperature and magnetic- field dependences of the specific heat, 
magnetization and susceptibility have been calculated by changing M — 1, 2, 3 and oo, 
results for M = oo corresponding to those of the conventional Boltzman-Gibbs statistics 
(BGS). It has been shown that the thermodynamical property of nanoclusters containing 
a small number of dimers is considerably different from that of macroscopic counterparts 
calculated within the BGS The specific heat and susceptibility of spin dimers described 
by the Heisenberg model have been discussed also by employing the NES. 
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1 Introduction 

1.1 A brief review of nanomagnetism 

In the last decade, there has been a considerable interest in atomic engineering, which 
makes it possible to create small-scale materials with the use of various methods (for 
reviews, see Refs. [1, 2, 3]). Small-scale magnetic systems ranging from grains (micros), 
nanosystems, molecular magnets and atomic clusters, display a variety of interesting prop- 
erties. Magnetic nanosystems consist of small clusters of magnetic ions embedded within 
nonmagnetic ligands or on nonmagnetic substrates. Nanomagnetism shows interesting 
properties different from bulk magnetism. Nanoclusters consisting of transition metals 
such as Fe N (iV=15-650) [4], Co^v (iV=20-200) [5], and Ni N (iV=5-740) [6] have been 
synthesized by laser vaporization and their magnetic properties have been measured, 
where N denotes the number of atoms per cluster. Magnitudes of magnetic moments 
per atom are increased with reducing N [6]. It is shown that magnetic moments in 
Co monatomic chains constructed on Pt substrates are larger than those in monolayer 
Co and bulk Co [7]. Recently Au nanoparticles with average diameter of 1.9 nm (in- 
cluding 212 atoms), which are protected by polyallyl amine hydrochloride (PAAHC), 
are reported to show ferromagnetism while bulk Au is diamagnetic [8]. This is similar 
to the case of gas-evaporated Pd fine particles with the average diameter of 11.5 nm 
which show the ferromagnetism whereas bulk Pd is paramagnetic [9]. The magnetic 
property of four-Ni molecular magnets with the tetrahedral structure (abbreviated as 
Ni4) in metallo-organic substance [Moi 2 O 30 (/i2 — OH) 10 H 2 {Ni(H 2 O3)}4] • 14H 2 has been 
studied [10]. Their temperature-dependent susceptibility and magnetization process have 
been analyzed by using the Heisenberg model with the antiferromagnetic exchange cou- 
plings between Ni atoms [10]. Similar analysis has been made for magnetic molecules 
of FeiV (N = 6, 8, 10 and 12) [11] [12], and V6 [13]. Extensive studies have been made 
for single molecule magnets of Mnl2 in [Mn 12 Oi 2 (CH 3 COO)i 6 (H 2 0) 4 ] [14] and Fe8 in 
[Feg(tanc)60 2 (OH) 12 ]Br9 ■ 9H 2 [15]. Both Mnl2 and Fe8 behave as large single spins 
with S = 10, and show quantum tunneling of magnetization and the square-root relax- 
ation, which are current topics in namomagnetism. Much attention has been recently paid 
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to single molecule magnets which are either dimers or behave effectively as dimers, due 
to their potential use as magnetic storage and quantum computing. The iron S = 5/2 
dimer (Fe2) in [Fe(OMe)(<i&m)2]2 [16] has a nonmagnetic, singlet ground state and its 
thermodynamical property has been analyzed with the use of the Heisenberg model [17]- 
[19]. Similar analysis has been made for transition-metal dimers of V2 [20], Cr2 [21], Co2 
[22], Ni2 [24] and Cu2 [25]. 

1.2 Non-extensive statistics 

As the size of systems becomes smaller, effects of fluctuations and contributions from 
surface play more important roles. There are currently three approaches to discussing 
nanothermodynamics for small-size systems: (1) a modification of the Boltzman-Gibbs 
statistics (BGS) adding subdivision energy [26], (2) non-equilibrium thermodynamics in- 
cluding work fluctuations [27], and (3) the non-extensive statistics (NES) generalizing 
the BGS as to take account of the non-extensive feature of such systems [28]-[33]. A 
comparison between these approaches have been made in Refs. [33] [34]. 

Before discussing the NES, let's recall the basic feature of the BGS for a system with 
internal energy E and entropy S, which is immersed in a large reservoir with energy E 
and entropy So- The temperature of the system T is the same as that of the reservoir T 
where T = SE/SS and To = SEo/SSq. If we consider the number of possible microscopic 
states of Q(Eo) in the reservoir, its entropy is given by So = kslnQ^Eo) where ks denotes 
the Boltzman constant. The probability of finding the system with the energy E is given 
by p(E) = tt(E - E)/Vt{E) ~ exp(-E/k B T) with E < E . When the physical quantity 
Q of a system containing N particles is expressed by Q oc TV 7 , they are classified into two 
groups in the BGS: intensive (7 = 0) and extensive ones (7 = 1). The temperature and 
energy are typical intensive and extensive quantities, respectively. This is not the case in 
the NES, as will be shown below. 

When a small-scale nanosystem is immersed in a reservoir, the temperature of the 
nanosystem is expected to fluctuate around the temperature of the reservoir T because 
of the smallness of the nanosystem and its quasi-thermodynamical equilibrium states. 
Then the BGS distribution mentioned above has to be averaged over the fluctuating 
temperature. This idea has been expressed by [31] [32] [33] 
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A) = rV= r d(3f{(5)(5 = E{(5), (4) 
2 £(/J 2 ) - £(,3) 2 



AT £(/3) 2 
where exp (x) denotes the g-exponential function defined by 



(5) 



exp g (a;) = [1 + (1 — q)x] , for 1 + (1 — q)x > 

= 0. otherwise (6) 

In Eqs. (l)-(6), q expresses the entropic index, f B (f3) the Y (or x 2 ) distribution function 
of the order N, E(Q) the expectation value of Q averaged over f(/3), (3o the average of 
the fluctuating f3 and 2/N its variances. The V distribution of the order N is emerging 
from the sum of squares of N Gaussian random variables. In deriving Eqs. (l)-(5), we 
have assumed that N particles are confined within a small volume of L 3 (L < £) where 
the variable (5 uniformly fluctuates, £ standing for the coherence length [32]. 

The important consequence of the NES is that energy and entropy are not proportional 
to N in nanosystems. The non-extensivity of the entropy was first demonstrated by 
Tsallis, who proposed the generalized entropy given by [28] 

s q = k B fil" 1 ) = - k B E p? MpO> (t) 

where Pi [= p(ei)] denotes the probability distribution for the energy in the system and 
ln ? (o;) [= (x l ~ q — 1)/(1 — q)\ the g-logarithmic function, the inverse of the g-exponential 
function defined by Eq. (6). It is noted that in the limit of q — 1, Eq. (7) reduces to the 
entropy of BGS, Sbg, given by 

S x = S B g = -k B ^Pi\v-Pi. (8) 

i 

The non-extensivity in the Tsallis entropy is satisfied as follows. Suppose that the total 
system containing 2N particles is divided into two independent subsystems, each of which 
contains N particles, with the probability distributions, p^p and pf\ The total system 
is described by the factorized probability distribution pij = p^ p^p . The entropy for the 
total system S(2N) is given by [28] 

S(2N) = S(N) + S(N) + O (1) , (9) 

where S(N) stands for the entropy of the iV-particle subsystem, the index q given by Eq. 
(2) being employed. Similarly the energy of the total system is expressed by 

E(2N) = E(N) + E(N) + O (1) , (10) 

The difference of E(2N) — 2E(N) is attributed to the surface contribution. This implies 
that the index 7 in Q oc iV 7 is neither nor 1 for Q = S and E in nanosystems within 
the NES. 
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The functional form of the probability distribution p(E) expressed by Eq. (1) was 
originally derived by the maximum-entropy method [28] [29]. The probability of [= p(e,)] 
for the eigenvalue q in the NES is determined by imposing the variational condition to 
the entropy given by Eq. (7) with the two constraints [29]: 

Eft = x > ( n ) 

i 

^? = E q . (12) 

The maximum-entropy method leads to the probability distribution pi given by 

Pi oc exp [-fa (a - E q )] , (13) 



with 



A) = -, (14) 



c, 



q 



EpI (15) 



where /3 denotes the Lagrange multiplier relevant to the constraint given by Eq. (12). It 
has been shown that the physical temperature T of the nanosystem is given by 



T — 7~5> (AMP) (16) 

KbP 

In the limit of q — 1, we get expjx] = e x , c q — 1 and given by Eqs. (13)-(16) reduces 
to the results obtained in the BGS, related discussions being given in Sec. 4. 

In previous papers [36, 37, 38], I have applied the NES to the Hubbard model, which 
is one of the most important models in solid-state physics (for a recent review, see Ref. 
[39]). The Hubbard model consists of the tight-binding term expressing electron hoppings 
and the short-range interaction between two electrons with opposite spins. The Hubbard 
model provides us with good qualitative description for many interesting phenomena such 
as magnetism, electron correlation, and superconductivity. In particular, the Hubbard 
model has been widely employed for a study on transition-metal magnetism. In the limit 
of strong interaction (U/t <C 1), the Hubbard model with the half-filled electron occu- 
pancy reduces to the Heisenberg or Ising model. The two-site Hubbard model has been 
adopted for a study on some charge-transfer salts like tetracyanoquinodimethan (TCNQ) 
with dimerized structures [40]- [42]. Their susceptibility and specific heat were analyzed by 
taking into account the interdimer hopping within the BGS. The NES calculations have 
been made for thermodynamical properties of canonical [36] [38] and grand-canonical en- 
sembles [37] of Hubbard dimers, each of which is described by the two-site Hubbard model. 
It has been shown that the temperature dependences of the specific heat and susceptibil- 
ity is significantly different from those calculated by the BGS when the entropic index q 
departs from unity for small iV [Eq. (2)], the NES in the limit of q — 1 reducing to the 
BGS. 
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The purpose of the present paper is to show (1) how thermodynamical property of a 
nanocluster containing a small number of Hubbard dimers is different from that of macro- 
scopic systems, and (2) how thermodynamical property of a given nanocluster is changed 
when M, the number of Hubbard dimers contained in it, is varied. The paper is orga- 
nized as follows. In Sec. 2, I apply the NES to nanoclustes, providing expressions for the 
energy, entropy, magnetization, specific heat and susceptibility. Numerical calculations 
of the temperature and magnetic-field dependences of thermodynamical quantities are 
reported for various M values. The final Sec. 4 is devoted to discussions and conclusions. 
In the Appendix, the NES has been applied to a cluster containing spin dimers described 
by the Heisenberg model. 



2 Nonextensive thermodynamics of Hubbard dimers 
2.1 Energy and entropy 

I have adopted a system consisting of sparsely distributed N c nanoclusters, each of which 
contains independent M dimers. It has been assumed that the distance between nan- 
oclusters is larger than £, the coherence length of the fluctuating (3 field, and that the 
linear size of the clusters is smaller than £. Physical quantities such as the entropy and 
energy are extensive for N c , but not for M in general [32]. 
The Humiltonian of the cluster is given by 

M 

H = ]T Hf\ (17) 

2 2 

H ? = -* EfaLaar + alalia) + n tf n jl - A*b-B 2(nj T - n n ), 

(1,2 ei) (18) 

where denotes the two-site Hamiltonian for the £th dimer, n^ a = a^ a aj a , a^ a the 
annihilation operator of an electron with spin a on a site j (e £), t the hopping integral, 
U the intraatomic interaction, /ib the Bohr magneton, and B an applied magnetic field. 
In the case of the half-filled occupancy, in which the number of electrons is N e = 2, six 
eigenvalues of are given by 

e iE = 0, 2fi B B, -2fi B B, U, ^ + A, ^ - A, for i = 1 - 6, t = 1 - M (19) 

where A = y / f/ 2 /4 + At 2 [40] [42] . The number of eigenvalues of the total Hamiltonian H 
is6 M . 

First we employ the BGS, in which the canonical partition function for H is given by 
[40] [42] 

Z BG = Trexp(-/3#), (20) 
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6 6 

E • • 2 exp[-/3( eil + • • +e iM )}, 



(d) ]M 
BG\ 1 



7(d) 
J BG 



[z. 

1 + 2 cosh(2p fi B B) + e _/?t/ + 2 e - /3C//2 cosh(/?A), 



(21) 

(22) 
(23) 



where (3 = 1/ksT, Tr denotes the trace and Z^q the partition function for a single 
dimer. By using the standard method in the BGS, we can obtain various thermodynamical 
quantities of the system [40, 41, 42]. Because of a power expression given by Eq. (22), 
the energy and entropy are proportional to M: Ebg — ME^q and Sbg — MS^q where 
EgQ and S^q are for a single dimer. This is not the case in the NES as will be discussed 
below. 

Next we adopt the NES, where the entropy S q for the quantum system is defined by 
[28] [29] 

*=*»(W)- (24) 

Here p q stands for the generalized canonical density matrix, whose explicit form will be 
determined shortly [Eq. (27)]. We will impose the two constraints given by 



Tr (p q ) 
Tr (p«g) 
Tr (pi) 



1, 

< H > r . 



(25) 
(26) 



where the normalized formalism is adopted [29] . The variational condition for the entropy 
with the two constraints given by Eqs. (25) and (26) yields 



Pi 



X, 



-exp Q 



'H - E q ) 



with 



X q = Tr exp 



-i^\(H-E n 



c q = Tr (pi) = X]~\ 



(27) 



(28) 
(29) 



where exp q (x) is the g-exponential function given by Eq. (6) and (3 is a Lagrange multiplier 
given by 



1 dE n 



(30) 



The trace in Eq. (28) and (29) is performed over the 6 eigenvalues, for example, as 

6 6 / 



H=l «M=1 V 



. C„ ) 



E q ) 



- E q ) 



(31) 
(32) 
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where the following conventions are adopted: 

Q = Qi + ■ ■ + e i M ) (33) 

E = £•• E ■ (34) 

It is noted that in the limit of q — 1, Eq. (31) reduces to 

X x = Zrq exp^] = [zgl exp (35) 

For g 7^ 1, however, X g cannot be expressed as a power form because of the property of 
the g-exponential function: 

exp q (x + y) ^ exp q (x) exp q (y). (for q ^ 1) (36) 

It is necessary to point out that E q in Eq. (26) includes X q which is expressed by E q 
in Eq. (28). Then E q and X q have to be determined self-consistently by Eqs. (26)-(29) 
with the T — f3 relation given by Eq. (16) for a given temperature T. The calculation 
of thermodynamical quantities in the NES generally becomes more difficult than that in 
BGS. 



2.2 Specific heat 

The specific heat in the NES is given by [36] 



c = fd/3\ (dE q 



\*rj\dpj- (37) 

Because E q and X q are determined by Eqs. (26)- (29), we get simultaneous equations for 
dE q /d(3 and dX q /df3, given by 



with 



dE q 
~dj 

dXq 

d(3 



an 



a>2i 



'dE q 

JP 
'dEg 

dp 



Ol2 



+ a 2 2 



' dX q \ 

dp r 



bi, 



a n = qpX q q 2 J2 w i 9 ^ 

i 

a 12 = -X- x E q -Pq{q-\)X^Y. w ^~ X ^- E ^ 

i 

a 2 i = PXl 
a 2 2 = 0, 



h = 



Wi = exp q 



i 



- E q ) 



(38) 
(39) 

(40) 
(41) 

(42) 
(43) 
(44) 

(45) 
(46) 
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The specific heat is then given by 



C, = (I) ( * ) . ,47) 

\ dT / \ 1 - an - a 12 a 2 i J 



with 

d/3 / /3 2 



(48) 



dT \X\- q -P(l-q)Xt q {dX q ldP))' 
In the limit of q — > 1, Eqs. (38)-(46) yield the specific heat in the BGS, given by [37] 

Cbg = = k B (3\< e\ > x -<e t > 2 ), (49) 

where < • >i is defined by Eq. (26) with q — 1: 

< Qj > 1= X 1 - 1 ^exp[-/^(e l - ^)] = ^exp(-/3 Q ) Q f . (50) 

i i 

2.3 Magnetization 

The field-dependent magnetization m q in the NES is given by [36] 

a£ "+r'A 7 ^. (52) 



^ 9 9fi ' 

By using Eqs. (26)-(29), we get the simultaneous equations for dE q /dB and dX q /dB 
given by 

^ g dE q dX q 

-QB = ail W + ° 12 W + dl ' (53) 

dX q dE q dX q J 

= a 21 — ^ + a 22 — ^ + d 2 , (54) 



<9B <95 OB 



with 



d 1 = -X-^w^i + PqXf^w^eifH, (55) 

i i 

d 2 = fiX*- 1 ^*,?!*, (56) 

i 

where //j = —dei/dB, and a^- = 1,2) are given by Eqs. (40)-(43). From Eqs. 
(51)-(56), we obtain m q given by 

( -c l2 + ^X-^l- Cll ) \ 

m q = q - d 2 , (57) 

\ 1 - c u - c 12 c 21 J 

= X^J2^^=<^> q - (58) 
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In the limit of q — > 1, Eqs. (55) and (56) reduce to 



di = - (fx i ) 1 + P{e i iii) 1 , (59) 
d 2 = PX,^),, (60) 



where < • >i is given by Eq. (50). By using Eq. (58), we get 



m BG = (fi i ) 1 , (61) 
4// B sinh(2/?B) 

— rv 1 

^BG 

where Z BG and < • >i are given by Eqs. (20) and (50), respectively. 
2.4 Susceptibility 

The high-field susceptibility in the NES is given by 

X q {B) = ^. (63) 

The zero-field susceptibility Xq(B = 0) i s gi ven by [36] 

Xq = Xq (B = Q) = -E^+^ l X^Xf\ (64) 

where = d 2 E q /dB 2 \ B=Q and X® = d 2 X q /dB 2 | B=0 . With the use of Eqs. (26)- (29), 
we get simultaneous equations for E^ and X^ given by 

= a n Ef + a 12 Xf + A, (65) 

X® = a 2X Ef + a 22 Xf + / 2 , (66) 

with 

h = -2/?gX«- 2 ^>f-VL (67) 

h = (PqXp-^up-^l (68) 

i 

where = 1,2) are given by Eqs. (40)-(43). From Eqs. (64)-(68), we get 

x q = — =/3^r 2 E^ _ vnB=o. (69) 

°21 i 

In the limit of q — 1, Eq. (69) yields the susceptibility in BGS: 

Xbg = (3 < l4 \b=o>i, (70) 



k B T) 3 + e-W + 2e-? u / 2 cosh(/3A) ' 
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3 Calculated results 



3.1 Temperature dependence 

In order to study how thermodynamical quantities of a cluster containing M Hubbard 
dimers depend on M, I have made some NES calculations, assuming the M — q relation 
given by 

which is derived from Eq. (2) with M = 2N for dimers. Simultaneous equations for E q 
and X q given by Eqs. (26)-(29) have been solved by using the Newton-Raphson method 
with initial values of E\ and X 1 obtained from BGS (q = 1) corresponding to M = oo in 
Eq. (72). Calculated quantities are given per dimer. 

Figures 1(a), 1(b) and 1(c) show the temperature dependence of the specific heat C q 
for U/t — 0, 5 and 10, respectively, with various M values. The specific heat for M = oo 
shown by bold solid curves, expresses the result in BGS, and it has a peak at lower 
temperatures for the larger interaction, as previous BGS calculations showed [41]. Note 
that the horizontal scales of Fig. 1(c) are enlarged compared to those of Figs. 1(a) and 
1(b). The peak becomes broader for smaller M. 

The temperature dependence of the susceptibility x q f° r U/t = 0, 5 and 10 is plotted 
in Figs. 2(a), 2(b) and 2(c), respectively. The susceptibility for M = oo (BGS) shown by 
the bold solid curve, has a larger peak at lower temperatures for larger U [41]. Note that 
the horizontal and vertical scales of Fig. 2(c) are different from those of Figs. 2(a) and 
2(b). We note that for smaller M, the peak in Xq becomes broader, which is similar to 
the behavior of the specific heat shown in Figs. l(a)-l(c). 

When the M value is varied, maximum values of the specific heat (C*) and the sus- 
ceptibility (x* q ) are changed, and the temperatures (T£ and T*) where these maxima are 
realized, are also changed. Figure 3(a) depicts and T* for U/t = 5 as a function of 
1/M. It is shown that with increasing 1/M, T* is much increased than Tq. Similarly, 
the 1/M dependences of C* and x* for U/t = 5 are plotted in Fig. 3(b), which shows 
that maximum values of C q and Xq are decreased with decreasing M. This trend against 
1/M is due to the fact that a decrease in M (= 2iV) yields an increase in fluctuations of 
(3 fields, and then peaked structures of the specific heat and susceptibility realized in the 
BGS, are smeared out by (3 in Eq. (1). 



3.2 Magnetic-field dependence 

Next I discuss the magnetic-field dependence of physical quantities. Figure 4 shows the B 
dependence of the magnetization m q for U/t = 0,5 and 10 with M = 2 at k B T/t = 1. For 
U/t = 0, m q in the NES is smaller than that in the BGS at fisB/t < 1, but at fisB/t > 1 
the former becomes larger than the latter. In contrast, in cases of U/t = 5 and 10, m q 
in the NES is larger than that in the BGS for ngB/t > 0. In order to study the B 
dependence in more details, I show in Fig. 5 the B dependence of the six eigenvalues of 
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€i for U/t = 5 [Eq. (19)]. We note the crossing of the lowest eigenvalues of £3 and at 
the critical filed: 

flP U 

^Sc=y T g+* 2 - T , (73) 

leading to /i B B c /t = 0.351 for U/t = 5.0. At B < B c (B > B c ), e 6 (e 3 ) is the ground 
state. At B = B c the magnetization m q is rapidly increased as shown in Figs. 6(a) and 
6(b) for ksT/t = 1.0 and 0.1, respectively: the transition at lower temperatures is more 
evident than at higher temperatures. This level crossing also yields a peak in Xq [Figs. 
6(c) and 6(d)] and a dip in C q [Figs. 6(e) and 6(f). It is interesting that the peak of \q 
for M = 2 is more significant than that for M = 00 whereas the dip of C q for M = 2 
is broader than that for M = 00. When the temperature becomes higher, these peak 
structures become less evident as expected. Similar phenomenon in the field-dependent 
specific heat and susceptibility have been pointed out in the Heisenberg model within the 
BGS [43]. 

In the case of the quarter-filled occupancy (N e = 1), the eigenvalues are q = —t—fisB, 
—t+fisB, t—fisB, and t+fXsB for i = 1—4. Although the level crossing occurs between €2 
and €3 at hbB = t, it does not show any interesting behavior because the crossing occurs 
between the excited states. The case for the three-quarter-filled occupancy (N e = 3) is 
the same as that of the quarter-filled occupancy because of the electron-hole symmetry 
of the model. 

Figure 6(b) reminds us the quantum tunneling of magnetization observed in magnetic 
molecular clusters such as Mnl2 and Fe8 [14] , which originates from the level crossings of 
magnetic molecules when a magnetic field is applied [14]. 



4 Discussions and conclusions 

I have applied the NES to Hubbard dimers for a study of their thermodynamical proper- 
ties. The current NES is, however, still in its infancy, having following unsettled issues, 
(i) For relating the physical temperature T to the Lagrange multiplier f3, I have employed 
the T — f3 relation given by Eq. (16). There is an alternative proposal with the T — (3 
relation given by [29] 

T — t~~7T' (TMP) (74) 
KbP 

which is the same as in the BGS. At the moment, it has not been established which of the 
AMP and TMP methods given by Eqs. (16) and (74), respectively, is appropriate as the 
T — (3 relation in the current NES. It has been demonstrated that the negative specific 
heat of a classical gas model realized in the TMP method [44], is remedied in the AMP 
method [35] . Recent theoretical analyses also suggest that the AMP method is better than 
the TMP method [45] [46]. The TMP method yields an anomalously large Curie constant 
of the susceptibility in the free spin model [37] [47] and in the Hubbard model [36, 37]. In 
my previous papers [36]- [38], NES calculations have been made by using the TMP and 
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AMP methods. It has been shown that both methods yield qualitatively similar results 
although there are some quantitative difference between the two: the non-extensivity in 
the TMP method generally appears more significant than that in the AMP method, 
(ii) The N — q relation given by Eq. (2) was obtained in Eqs. (l)-(5) with the T distribu- 
tion f B ((3) given by Eq. (3). Alternatively, by using the large-deviation approximation, 
Touchette [48] has obtained the distribution function f T (/3), in place of f B (/3), given by 

fW = F{S) {—) " 2 exp {-W> ' (75) 

For N — > oo, both f B (j3) and f T {f3) distribution functions reduce to the delta-function 
densities, and for a large iV (> 100), both distribution functions lead to similar results. 
For a small iV (< 10), however, there is a clear difference between the two distribution 
functions (see Fig. 4 of Ref. [38]). It should be noted that f T cannot lead to the q- 
exponential function which plays a crucial role in the NES. For a large e, the T distribution 
f B in Eq. (1) yields the power form of w(e) ~ e~~ while f T substituted to Eq. (1) leads 
to the stretched exponential form of w(e) ~ e c ^. This issue of / versus f T is related to 
the superstatistics, which is currently studied with much interest [49]. 

To summarize, I have discussed thermodynamical properties of a nanocluster contain- 
ing M dimers, applying the NES to the Hubbard model. It has been demonstrated that 
the thermodynamical properties of a nanocluster with a small M calculated by the NES 
may be considerably different from those obtained by the BGS. It is interesting to com- 
pare our theoretical prediction with experimental results for samples containing a small 
number of transition-metal dimers. Unfortunately samples with such a small number of 
dimers have not been reported: samples having been so far synthesized include macro- 
scopic numbers of dimers, to which the present analysis cannot be applied. I expect that 
it is possible to form a dimer assembly by STM manipulation of individual atoms [50]. 
Scanning probes may be used also as dipping pens to write small dimerized structures [51]. 
Theoretical and experimental studies on nanoclusters with changing M could clarify a link 
between the behavior of the low-dimensional infinite systems and finite-size nanoscale sys- 
tems. I hope that the unsettle issues (i) and (ii) in the current NES mentioned above 
are expected to be resolved by future experiments on nanosystems with changing their 
sizes. It would be interesting to adopt quantum-master-equation and quantum-Langevin- 
equation approaches, and/or to perform large-scale molecular-dynamical simulations, for 
nanoclusters described by the Hubbard model. 
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Appendix: NES for Heisenberg dimers 

I have considered a cluster containing M spin dimers (called Heisenberg dimers) de- 
scribed by the Heisenberg model (s = 1/2) given by 



H 



M 

- Jsi • s 2 - gfi B B(s lz + s 2z ), 



[1,2 e 



(76) 
(77) 



where J stands for the exchange interaction, g (=2) the g-factor, fi B the Bohr magneton, 
and B an applied magnetic field. Four eigenvalues of are given by 



3J 



- g/j, B Bmi 



with mi = 1, 0, — 1 for i = 1,2, 3, 
with rri4 = for % = 4, 



In the BGS the canonical partition function is given by [17]- [19] 



J BG 



J BG 



— [Zbg^ 



7{d) iM 
'B 

exp 



^ J [1 + 2cosh(g fx B (3B)} + exp 



(78) 

(79) 
(80) 



with which thermodynamical quantities are easily calculated. The susceptibility is, for 
example, given by 

Xbg = Mxbg, 



(d) 
Xbg 



/4 



(81) 

k B T\3 + exp(-J/k B T)) ' (82) 

The calculation of thermodynamical quantities in the NES for the Heisenberg model 
goes parallel to that discussed in Sec. 2 if we employ eigenvalues given by Eq. (78). For 
example, by using Eq. (69), we get the susceptibility for the Heisenberg model, given by 



2 2 (qP\ 1 V- 

X q = 9 M — IT 2^ 

\ Q / ^1 i 

In the case of M — 1 (a single dimer), we get 



2o-l 2 



cxp ( . 



p' 



J 



- \ -v + E, 



2q-l 



with 



x„ 



E n 



3 exp. 



P' 



J 



- ^ + E, 



< 



-3J 



exp. 



P 



+ exp 
J 



P' 



3J F 
T 4 



(83) 



(84) 



(85) 



+ E C 



+ (t)( 6XP « 



I) f^L 

cJU 
14 



-E, 



(86) 



In the limit of q = 1, Eq. (84) reduces to Xbg gi ven by Eq. (82). 

The Curie constant Y q defined by Xq — (f 1 B/^B)(Xq/T) for T ^> J is given by 

T q = 2Mq, (AMP) (87) 

= 2MqA M(q - 1) . (TMP) (88) 

Equations (87) and (88) are derived with the use of the T — f3 relation given by Eqs. (16) 
and (74), respectively. These are consistent with results obtained for Hubbard dimes [37]. 

Figures 7(a) and 7(b) show the temperature dependence of the specific heat C q and 
susceptibility Xq °f Heisenberg dimers calculated with the use of Eq. (83) for M — 1, 2, 
3 and oo (M = oo corresponding to the BGS with q = 1.0). We note that the results 
of Heisenberg dimers are quite similar to those of the Hubbard dimer for U/t = 5 and 
10 shown in Figs. 2(b), 2(c), 2(e) and 2(f). This is not surprising because the Hubbard 
model with the half-filled electron occupancy in the strong-coupling limit reduces to the 
Heisenberg model. 
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Figure 1: The temperature dependence of the specific heat C q per dimerior (a) U/t = 0, 
(b) 5 and (c) 10, calculated for M = 1 (solid curves), 2 (chain curves), 3 (dashed curves) 
and oo (bold solid curves), results for M = oo denoting those in the BGS. 



Figure 2: The temperature dependence of the susceptibility Xq P er dimerior (a) U/t = 0, 
(b) 5 and (c) 10, calculated for M = 1 (solid curves), 2 (chain curves), 3 (dashed curves) 
and oo (bold solid curves), results for M = oo denoting those in the BGS. 



Figure 3: (a) 1/M dependence of the temperatures of Tq (circles) and T* (squares) where 
C q and Xq have the maximum values, respectively, (b) 1/M dependence of the maximum 
values of C* (circles) and x* q (squares) (U/t = 5) 



Figure 4: The magnetic-filed dependence of the magnetization m q for (a) U/t — 0, (b) 5, 
and (c) 10 with k B T/t = 1 for M = 2 (solid curves) and oo (dashed curves). 



Figure 5: The magnetic-filed dependence of the eigenvalues (i — 1 — 6) for U/t — 5, B c 
denoting the critical field where a level crossing between and e$ occurs. 



Figure 6: The magnetic-filed dependence of (a) the magnetization m q for k B T/t = 1.0 
and (b) ksT/t = 0.1, (c) the susceptibility Xq f° r ksT/t = 1.0 and (d) k B T /t = 0.1, and 
(e) the specific heat C q for k B T/t = 1.0 and (f) k B T/t = 0.1 with U/t=5, calculated for 
M = 2 (solid curves) and oo (dashed curves). 



Figure 7: The temperature dependence of (a) the specific heat and (b) susceptibility of 
Heisenberg dimers for various M: M — 1 (bold solid curves), 2 (chain curves), 3 (dashed 
curves), and oo (solid curves). 
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